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Detection of LS I +61 303 in a low VHE 7-ray emission state with the MAGIC telescopes 
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Abstract: The gamma-ray binary system LS I +61 303 was studied in great detail in VHE gamma-rays in the last years 
by the MAGIC telescope. The VHE emission of the system exhibited a prominent periodic outburst in the orbital phases 
0.6-0.7 between September 2005 to lanuary 2008. In Fall 2008 the Fermi collaboration reported as well periodic emission 
in the MeV to GeV energy range, but with a shifted outburst in the phases 0.35-0.45. MAGIC observed again LS 1+61 303 
in 2009 with the twice more sensitive stereo mode to allow for detailed correlation studies between the VHE gamma-ray 
and Fermi/LAT energy band. Here we present our new results, which show a significant reduction in the VHE gamma- 
ray flux in the phase of the periodic outburst by almost one order of magnitude compared to our previous measurements. 
Furthermore, the 0.1-phase averaged light curve shows no significant outburst, but a rather constant flux. Here we will 
discuss the implications of our results for future gamma-ray studies of LS I +61 303. 
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1 Introduction 

The LS I +61 303 system consists of a Be star and a com- 
pact object of still uncertain nature, either a neutron star 
or a black hole. Its orbital period, which is most precise 
measured in radio, is 26.4960 ± 0.0028 days |1|. Soft 
X-ray outbursts modulated with the same period as in the 
radio waveband were reported by [2| but the peak of the 
emission changes its phase over several years 0. Many 
orbital parameters of the system remain controversial (see 
G]|5] [4)) but observations indicate a highly eccentric orbit 
(e = 0.55 ± 0.05) with the periastron passage at orbital 
phase pcr = 0.275 [4|. These orbital parameters are im- 
portant for modelling the VHE emission of the system as 
shown in,e.g., Q or (8). 

In 2006 the MAGIC collaboration discovered variable 
VHE 7-ray emission from LS I +61 303 HD. A follow- 
ing extensive observational campaign in Fall 2006 found 
a period for the VHE emission of 26.6 ± 0.2 days |16|. 
The VHE 7-ray emission showed an outburst in the orbital 
phase interval 0.6-0.7 with no significant 7-ray emission 
detected during the rest of the orbit. In particular, no VHE 
7-ray signal was detected by MAGIC around the perias- 
tron passage of the system. The data from Fall 2006 also 
suggested a correlation between the X-ray and VHE 7-ray 
flux. An extensive multi-wavelength campaign conducted 
in 2007, including MAGIC, XMM-Newton and Swift, pro- 
vided strong evidence for the X-ray/VHE 7-ray flux cor- 
relation in the strictly simultaneous taken data 1 1 8 1 . In 
contrast, no correlation was found between the radio wave- 



length flux and the VHE 7-ray flux from the Fall 2006 cam- 
paign 03). 

The VHE emission of LS I +61 303 was confirmed by 
VERITAS [ 1 1 . However, in observations conducted by 
the VERITAS collaboration in Fall 2008 and early 2009, 
no VHE signal was detected. More recent VERITAS ob- 
servations in Fall 2009 (the same time period considered 
in the present paper) also yielded only upper limits for the 
VHE emission from LS I +61 303 Oil- Very recently the 
VERITAS collaboration reported a detection of the system 
with a significance of more than 5cr ( ll26irTT1l ) between or- 
bital phases 0.05 and 0.23. This places the detection at su- 
perior conjunction and 5.8 to 1.3 days before the periastron 
passage a phase range where no VHE 7-ray emission was 
previously detected. 

The binary system was observed in high energy (HE, 0.1 — 
lOOGeV) 7-rays by EGRET (23, but the large position un- 
certainty of the source and inconclusive variability studies 
of the emission|31 1, prevented its unambiguous identifica- 
tion. The spacial association with LS 1 +61 303 was only 
achieved following HE 7-ray observations by AGILE [28 ]. 
Most recently, Fermi/LAT found that the HE 7-rays are pe- 
riodically modulated in very good agreement with the (ra- 
dio) orbital period [9] and firmly established LS I +61 303 
as the counterpart of the HE emission. The HE outburst 
was not, however, observed at the same phases as the VHE 
outburst but occurred just after the periastron passage be- 
tween phases 0.3-0.45. This difference in phase may indi- 
cate that different processes are responsible for the HE and 
VHE 7-ray emission. On the other hand, the same process 



T. JOGLER ET AL. LS I +61 303 IN A LOW VHE 7-RAY EMISSION 



might produce both emissions if the GeV 7-rays are pro- 
duced by inverse Compton (IC) pairs cascading developing 
in the radiation field of the star. For more details on such 
a scenario see [20]. Another possibility is that the shift in 
the peak emission could be caused by a different location 
of the 7-ray production site in the system [32|. We note 
that no simultaneous VHE observations are available for 
the same time (Aug 2009 to Jan 2009) of the first reported 
Fermi observations |9|. An unambiguous interpretation of 
the non simultaneous SED from MeV to TeV energies of 
LS I +61 303 is not possible because the system might have 
changed its VHE emission in the meantime. 

Two principal scenarios have been proposed to explain the 
nonthermal emission from LS I +61 303: an accretion pow- 
ered microquasar (e.g. |30l [T9j |25] |22l) and a rotational 
powered compact pulsar wind (e.g. Il24l [7] [331 ). An al- 
ternative model assumes that the compact object is an ac- 
creting magnetar and that the 7-rays are produced along 
the accretion flow onto the magnetar OTTl . Although high 
resolution radio measurements ll23l show an extended and 
variable structure emerging from the system, it is not clear 
if they are produced by collimated plasma outflows (jets) or 
by a pulsar wind interacting with that of the Be star. Neither 
of the two proposed scenarios could be validated by accre- 
tion disk features, e.g. a thermal component in the X-ray 
spectrum, or the presence of pulsed emission at any wave- 
length. Thus the engine behind the VHE emission remains 
an open question. 

Here we present new observations of LS 1 +61 303 con- 
ducted with the MAGIC stereo system. This has twice 
the sensitivity of the previous MAGIC campaigns, and re- 
sults in a significant detection of the binary system during 
a newly identified low emission state. 

2 Observations 

The observations were performed between 2009 Oct 15 and 
2010 Jan 22 using the MAGIC telescopes on the Canary is- 
land of La Palma (28.75°N, 17.86°W, 2225 m a.s.l.), from 
where LS I +61 303 is observable at zenith distances above 
32°. The MAGIC stereo system consists of two imaging air 
Cherenkov telescopes, each with a 17 m diameter mirror. 
The observations were carried out in stereo mode, mean- 
ing that only shower images which trigger simultaneously 
both telescopes are recorded. The stereoscopic observa- 
tions provides a 5er signal above 300 GeV from a source 
which exhibits 0.8% of the Crab Nebula flux in 50 hours 
observation time, a factor of two more sensitive than our 
single telescope campaign on LS 1 +61 303 in 2007. Fur- 
ther details on the design and performance of the MAGIC 
stereo system can be found in IfTTll . 

The LS I +61 303 data set spans four orbits of the system, 
with two of these orbits observed for only one and three 
nights, respectively. The data taken in 2009 Oct and 2009 
Nov were restricted to moonless nights. The data sample 
the orbital phases 0.55 to 0.975 for 2009 Oct , and 0.575 



to 1.025 for 2009 Nov , the last night of which is in the 
next orbital cycle. The data recorded in 2010 Jan cover the 
phases 0.22 to 0.32 and were obtained during moonlight 
conditions. All data were taken at zenith angles between 
32° and 48°. After pre-selection of good quality data a 
total of 48.4 hours of data remained for the analysis. The 
observation strategy aimed to cover consecutive nights with 
at least three hours of observation in each individual night. 
Due to adverse observation conditions such as bad weather 
or too high zenith angle, the data set does not have uniform 
coverage during the orbital phases and some nights have 
shorter observation times than the planned three hours. 

3 Data Analysis 

The data analysis was performed with the standard MAGIC 
reconstruction software. The recorded shower images were 
calibrated, cleaned and used to calculate image parame- 
ters individually for each telescope. The energy of each 
event was then estimated using look up tables generated by 
Monte Carlo (MC) simulated 7-ray events. Events that trig- 
gered both telescopes simultaneously were then selected 
(so called stereo events^ and further parameters, e.g. the 
height of the shower maximum and the impact param- 
eter from each telescope, were calculated. The gamma 
hadron classifications and reconstructions of the incoming 
direction of the primary shower particles were then per- 
formed using the Random Forest (RF) method |14|. Fi- 
nally, the signal selection used cuts in the hadronness (cal- 
culated by the RF) and the squared angular distance be- 
tween the shower pointing direction and the source position 
(9 2 ). The energy dependent cut values were determined by 
optimizing them on a sample of events recorded from the 
Crab Nebula under the same zenith angle range and similar 
epochs than the LS I +61 303 data. For the energy spec- 
trum and flux, the effective detector area was estimated by 
applying the same cuts used on the data sample to a sample 
of MC simulated 7-rays. Finally, the energy spectrum was 
unfolded, accounting for the energy resolution and possible 
energy reconstruction bias |[T3l . 

4 Results 

The integral data set of 48.8 hours shows a 6.3<t detection 
of VHE 7-ray emission above 400 GeV from LS I +61 303 
(see Fig.[TJ. The integrated flux above 300 GeV is 

F(E > 300GcV) = (1.4±0.3 st ±0.4 sy )xl0" 12 cm- 2 s- 1 . 

(1) 

corresponding to about 1.25% of the Crab Nebula flux in 
the same energy range. 

1 . This step is only needed for the 2009 October data where no 
hardware stereo trigger was yet available. 
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Figure 1 : The squared angular distance between the point- 
ing direction of the shower and the source position (6> 2 -plot) 
for the position of LS 1 +61 303 (points) and the simul- 
taneous determined background regions (grey shaded his- 
togram) for the total MAGIC data set. 

4.1 Light curve 

We derived a nightly light curve above an energy of 
300 GeV that is shown in Fig. [2] A constant flux fit to the 
light curve yields a x 2 / dof = 42.15/19 (p = 1.5 x 10" 3 ) 
and hence is unlikely. Thus, as in previous observa- 
tions, the emission is variable and reaches a maximum 
flux around orbital phase 0.62 of F(E > 300GeV) = 
(6.1 ± 1.4 st ± l-8 sy ) x 10~ 12 cm _2 s _1 corresponding to 
5.4% of the Crab Nebula flux. This is a much lower peak 
emission than detected in our previous campaigns at the 
same orbital phases and sampled with very similar cadence. 
For a more quantitative comparison of the 2009 emission 
level with the previous MAGIC observations, we included 
the phase averaged 2007 light curve in Fig. [2] 

It is clear that the averaged emission level is lower than the 
data from 2005 to 2007. Not only had the flux changed but 
a VHE 7-ray excess was also observed at phases other than 
those of the periodic outburst between 0.6-0.7. The highest 
flux is, however, again detected in an outburst during the 
interval 0.6-0.7, and the measurements in the orbital cycle 
of 2009 Nov show the same burst profile as in previous 
observations but with a reduced flux level. Whether this 
outburst is still a periodic property of the VHE emission for 
LS I +61 303 cannot be determined with the small sample 
of orbital cycles observed in this campaign; it is noteworthy 
that the outburst was not detected during the orbit observed 
in 2009 Oct. 

The rather low mean values, even in the phase bin 0.6-0.7, 
of the individual night peak emission indicates that most 
of the emission of the system is contributed by only few 
nights. This is the case at least for the phase bins 0.6-0.7 
and 0.9-1.0, where the mean flux shows the highest signif- 
icance. On the other hand, there are similar fluxes, mea- 
sured in phase intervals 0.2-0.3 and 0.5-0.6, although with 
slightly lower significances, without any obvious individ- 
ual high flux nights. 
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Figure 2: VHE (E > 300GeV) 7-ray flux of LS I +61 303 
as a function of the orbital phase for the four observed or- 
bital cycles (four upper panels) and averaged for the en- 
tire observation time (lowermost panel, black points). The 
starting MjD for each orbital cycle is given in each panel. 
In the lower most panel the previous published [18] aver- 
aged fluxes per phasebin are shown as hollow triangles too. 
Vertical error bars show lcr statistical errors. 

4.2 Spectrum 

The emission level of LS I +61 303 was too low during 
most phases to obtain statistical significant phase depen- 
dent spectra. The total signal, however, was sufficient to 
form a phase integrated spectrum for the source. The spec- 
trum is well described by a simple power law 



dF 
dE 



(2.3±0.6 st ±0.7 s 



10 



-13 



E 



TeV cm 2 s 



ITeV 



-2.5±0.5 3t ±0.2 a - 



(2) 



with a reduced \ 2 = 0.24/2. The spectral fit parameters 
match those previously reported by MAGIC 02] [T6j [18) 
showing no evidence for long term spectral variability de- 
spite very different fluxes during these different campaigns. 

5 Conclusion 

The binary system LS 1 +61 303 was detected emitting 
VHE 7-ray s in 2009 at a level a factor 10 lower than pre- 
viously observed in the phase interval 0.6-0.7. The pre- 
viously observed orbital modulated outburst in this inter- 
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val was not detected during the first observed orbital cycle, 
whereas it was observed in the second. From this data set 
alone it was not possible to determine whether the outburst 
is still a truly periodic feature of the light curve. Although 
we could not significantly detect emission at individual or- 
bital phases in the here presented observations because of 
the weakness of the source, it appears that several phase 
intervals dominate in the integral signal. This is strong ev- 
idence for a new behavior in the VHE 7-ray emission of 
LS I +61 303. In previous observations conducted with 
MAGIC, the 0.6-0.7 interval dominated the total flux. 

Furthermore it is evident that the flux during the phase in- 
terval 0.6-0.7 is considerable reduced compared to the pre- 
vious campaigns, and on a similar level as in other phase 
intervals (e.g. 0.9-1.0). This suggests that a change in the 
VHE 7-ray emission of LS I +61 303 has occurred. On 
the other hand, there was no statistically significant change 
in the spectrum of the orbit-integrated flux in 2009 com- 
pared to the earlier results, suggesting that the same pro- 
cesses continue to produce VHE gamma-rays, but that ei- 
ther fewer are produced or they are more absorbed. 

If enhanced absorption decreases the observed VHE flux, 
the photon energy will be redistributed to lower energies 
and thus might be visible as a flux increase at lower ener- 
gies. The produced VHE 7-rays would thus need to prop- 
agate through a circum-source environment denser in pho- 
tons and we would expect to detect a cut-off or an absorp- 
tion feature in our spectrum, neither of which is found. De- 
spite the richness of observations of the system, its VHE 
emission is still not understood. 

The VERITAS observations in the same period as the one 
considered here did not detect VHE 7-ray emission from 
the system ifTTl . Our measurements are not, however, in 
contradiction to those of VERITAS. Our longer integra- 
tion combined with a denser sampling of two orbital cycle 
yielded a fainter detection threshold than the expected VHE 
7-ray signal from LS 1 +61 303 from previous campaigns 
expected. Thus it is evident that a frequent sampling with 
long individual integrations is required not to miss weak 
emission from binary systems. 

This is the first VHE 7-ray detection of LS I +61 303 
in the era of the Fermi satellite. The faint emission at 
VHE 7-rays does not yet permit night by night correla- 
tion studies but do show that the emission in LS I +61 303 
has changed on longer timescale since 2007. More sensi- 
tive and even deeper VHE 7-ray observations should yield 
shorter timescale correlation studies. 
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